The Evolution of the Kinematics of Nebular Shells in 
Planetary Nebulae in the Milky Way Bulge^ 



Michael G. Richer, Jose Alberto Lopez, Maria Teresa Garci'a-Di'az, 
David M. Clark, Margarita Pereyra 
OAN, Instituto de Astronomia, Universidad Nacional Autonoma de Mexico, 
P.O. Box 439027, San Diego, CA 92143 
{richer, jal, tere, dmclark, mally}(§astrosen.unam.inx 

and 

Enrique Diaz-Mendez 
Dept. of Physics and Astronomy, Texas Christian University, Fort Worth, Texas, U.S.A. 76109 

e . d . mendezQt cu.edu 



ABSTRACT 

We study the hue widths in the [O III]A5007 and Ha hues for two groups of planetary nebu- 
lae in the Milky Way bulge based upon spectroscopy obtained at the Observatorio Astronomico 
Nacional in the Sierra San Pedro Martir (OAN-SPM) using the Manchester Echelle Spectro- 
graph. The first sample includes objects early in their evolution, having high H/3 luminosities, 
but [O III] A5007/H/3 < 3. The second sample comprises objects late in their evolution, with 
He II A4686/H/3 > 0.5. These planetary nebulae represent e volutionary phases preceeding and 
following those of the objects studied bv [Richer et al. ] (l20Q8h . Our sample of planetary nebulae 
with weak [O III]A5007 has a line width distribution similar to that of the expansion velocities of 
the envelopes of AGB star s, and shifted to sys tematically lower values as compared to the less 



evolved objects studied bv [Richer et al. ] (l2008l) . The sample with strong He II A4 686 has a line 
width distribution indistinguishable from that of the more evolved objects from iRicher et al. I 
( 20081 ). but a distribution in angular si ze that is systematic ally larger and so they are clearly 
more evolved. These data and those of [Richer et al. I ( 20081 ) form a homogeneous sample from 
a single Galactic population of planetary nebulae, from the earliest evolutionary stages until the 
cessation of nuclear burning in the central star. They confirm the long-standing predictions of 
hydrodynamical models of planetary nebulae, where the kinematics of the nebular shell are driven 
by the evolution of the central star. 

ISM: planetary nebulae (general)- 



Subject headings. 
Galaxy: bulge 

Introduction 



Hydrodynamical models of planetary nebu- 
lae have long predicted a particular kinematic 
evolution for the nebular shells, driven pri- 



-ISM: kinematics and dynamics — stars: evolution- 



'^The observations reported herein were acquired at the 
Observatorio Astronomico Nacional in the Sierra San Pedro 
Martir (OAN-SPM), B. C, Mexico. 



(e.g., 


Kwok et al. 1978: Kahn & West 19851 


Schmidt- VoiEt & KoDoen .,1987a, b.:,Breitschwerdt & Kahn 


1990; 


Kahn & Breitschwerdt Ill990; Marten & Schonberner 


1991: 


Mellcma "l994':'ViUaver et al. l2002l Perinotto et al. 


.2004; .Schonberner et al. .,2007 


) . Initially, the cen- 



tral stars are cool and their winds relatively slow. 
This wind interacts with the wind that the precur- 
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sor asymptotic giant branch (AGB ) star emitte d 



in a momentum-conserving mode (jKwok 1 119821 ) . 
However, the central star's temperature and wind 
velocity increase rapidly, with the consequences 
that an ionization front is driven through the 
AGB envelope and the interaction between the 
two winds switches to an energy driven mode, 
and a hot bubble is created behind the shocked 
wind. The ionization front first accelerates the 
AGB envelope, now seen as the rim of the plane- 
tary nebula. In time, once the internal pressure of 
the hot bubble exceeds that of the nebular shell, 
it further accelerates the nebular shell. Theoreti- 
cally, the latest phases of evolution are less clear, 
though the central star will cease nuclear energy 
generation, fade rapidly, and emit an ever-weaker 
wind, in principle allowing the inner part of the 
nebul ar envelope to backfill into the central cavity 



(e.g., Garci'a-Segura et al. 2006[ l. 



Although the many extant observations of the 
kinematics of planetary nebulae all show expan- 
sion of the nebular shell, there are few systematic 
observations of how these shells acquir e their mo- 
tion and how it evolves with time. Dopita et al 



( 1985 . 19881) were the first to provide observational 
support for the early acceleration of the nebu- 
lar shell from studies of planetary nebulae in the 
Magellanic Clouds. Studies of Milky Way plane- 
tary nebulae prov ided much less convincing results 
(e.g., IChu et al. iri984,;.Heap . 1.993 ; ^Medi na et al 



2ooe 



^. RecentlvT lRicher et al. I (|20C)8t) demon- 



strated the acceleration of nebular shells in bright 
planetary nebulae in the Milky Way bulge during 
the early evolution of their central stars and were 
able to associate the acceleration seen in different 
evolutionary stages with the phases of acceleration 
expected from theoretical models. 

Here, we underta ke a study that complements 



Richer et al. 



( 20081 ). selecting objects earlier and 
later in their evolution than they did. The addi- 
tion of these objects allows us to study the evolu- 
tion of the kinematics of the nebular shell from the 
earliest stages of the planetary nebula phase un- 
til the cessation of nuclear burning in the central 
stars. In section [21 we present our new data and 
their analysis. In section [Sj our results and their 
implications are outlined, the principal ones being 
the similarity of the line widths in Ha and [O III] 
A5007, that our sample of least evolved objects 
has a line width distribution shifted to the lowest 



values while the sample of most evolved objects 
has a size distribution with the largest sizes, and 
that the evolutionary state correlates with the H/3 
luminosity. In section|3]we summarize our conclu- 
sions. 

2. Observations and Analysis 

2.1. The Planetary Nebula Sample 

Table[T]lists our two new samples of Bulge plan- 
etary nebulae. There are 24 objects in the sample 
with low [O III]A5007/H/3 ratios, all drawn from 



1987: Webster' 


'1988: 'Cuisinier et al. Ill996'. '2000|; 


RataE et al. .,1997.;.Escudero & Costa 12001: .Escudero et al. 


2004; 


Exter et al. 


2004 Gornv et al. 


20041). 



There are 21 objects in the sample with strong 
He II 4686, only some of which have extensive 
spectroscopy (see the previous references). Many 
of these objects were selected o n the basis o f 
the H e II 4686 line intensity from iTvlenda et al. 



(jl994l ) . Figure [T] explains the logic of the selection 
criteria for the two samples. 

The sample with weak [O III] A5007 com- 



prises objects entirely excluded from lRicher et aL 
(l2008l) and its selection were (i) a position 
within 10° of the galactic centre, (ii) a large ob- 
served, reddening-corrected H/3 flux, nominally 
log /(H/3) > — 12.0dex, and (iii) an intensity ratio 
[O III]A5007/ H/3 < 3. Since even t he less evolved 
objects from [Richer et al. ] (Sooi) show the ef- 
fects of acceleration of the nebular shell due to 
the passage of the ionization front, our last crite- 
rion should select objects that are significantly less 
evolved and thus probe the nebular kinematics at 
a stage that most closely reflects the kinematics of 
the undisturbed AGB envelope. The weak [O III] 
A5007 sample should include bright, young plane- 
tary nebulae whose central stars are still relatively 
cool. All of the central stars should be on the hor- 
izontal part of their post- AGB evolutionary track. 

The selection criteria for the "strong He II A4686" 
sample were (i) a position within 10° of the galac- 
tic centre and (ii) an intensity ratio He II A4686/H/3 > 
0.5. There was no restriction on the H/3 fiux. 
The second criterion should select a sample of 
objects biased to more advanced evolutionary 
ph ases than the more evolved objects included 



[Richer et al. ] (|2008[ ). since their requirement 



for their more evolved objects was equivalent to a 
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ratio He II A4686/H/3 > 0.09. This sample should 
include planetary nebulae whose central stars are 
at or slightly before their maximum temperature 
or fading towards the white dwarf cooling track. 

2.2. Observations, Data Reduction, and 
Derived Parameters 

We acquired the observations reported here and 
measured the derived properties in a very similar 
fashio n to our previous observations ([Richer et al 
2008 ). More details of our a nalysis may be found 



in [Richer et al. I (|2008l boOOl ) . 

We obtained high resolution spectra with 
the Manchester echel l e spec trometer (MES-SPM; 



Meaburn et al. 111984 120031 ) on 2008 June 1-2 and 
10-17 and 2008 September 1-11 at the Obser- 
vatorio Astronomico Nacional in the Sierra San 
Pedro Martir, Baja California, Mexico (OAN- 
SPM). A 150 ^m wide slit (1'.'9 wide on the sky, 
5' long) yielded a spectral resolution equivalent to 
llkm/s (2.6 pix FWHM) and a spatial sampling 
of O'.'6/pixel when coupled to a SITe 1024 x 1024 
CCD with 24 /^m pixels binned 2x2. We used a 
ThAr lamp for the wavelength calibration, which 
typically yielded an internal precision better than 
±1.0km/s. 

We usually obtained a single deep spectrum in 
each of the [O III] A5007 and Ha filters of, at most, 
30 minutes duration. When possible, the expo- 
sure time for the Ha and [O III] A5007 spectra 
were chosen to achieve similar signal levels. The 
slit was oriented north-south and, except for M 2- 
38, centered on the object. All of the planetary 
nebulae we observed are resolved (see Tabled]). 

The spectra were reduced using the twod- 
spec and specred packages of the Image Reduc- 
tion and Analysis Facilitji (IRAF), following 
Massev et aTl (|l992l Appendix B). We edited the 
spectra to remove cosmic rays and subtracted a 
nightly mean bias image. We rectified the object 
spectra and calibrated them in wavelength using 
the ThAr spectra. Finally, wavelength-calibrated, 
one-dimensional spectra were extracted for each 
object. No flux calibration was performed. 

The single exception to the foregoing was M 



^IRAF is distributed by the National Optical Astronomical 
Observatories, which is operated by the Associated Uni- 
versities for Research in Astronomy, Inc., under contract 
to the National Science Foundation. 



2-38, for which two 30 minute Ha spectra were 
obtained at different positions. The two spectra 
were reduced as described above and the final one- 
dimensional spectra were summed to produce the 
final Ha spectrum. 

We analyzed the one-dimensional spectra using 
a locally-implemen ted software package (INTENS; 
McCah et allll985l) . INTENS models the emission 
line with a sampled gaussian function and models 
the continuum as a straight line. For the strong 
He II A4686 sample, the Ha line was usually ac- 
companied by the He II A6560 line. In these cases, 
a fit was made simultaneously to both lines and 
the continuum, assuming that both lines had the 
same width. 

Table [T] presents the observed line widths 
(FWHM; full width at half the maximum inten- 
sity) for each object in both Ha and [O III] A5007. 
The formal uncertainties from fitting a Gaussian 
function with INTENS (one sigma; Table [T|) in- 
crease as the line width increases, a result of the 
line shape departing more from the Gauss ian form 
at larger line widths (jRicher et al. ] |2009l) . To ob- 



tain an idea of the real uncertainties, we measured 
the FWHM of the Ha line profiles directly (using 
implot/IRAF). For FWHM < 1.2A, there is no 
systematic difference between the line width and 
the Gaussian fit, though there is a dispersion of ap- 
proximately ±5% of the line width. For larger line 
widths, the Gaussian fit systematically underesti- 
mates the line width, with the difference reaching 
8 - 9% of the line width at FWHM - 2A. 

We derive the true, intrinsic profile {a true), re- 
sulting from the kinematics of the planetary neb- 
ula, by correcting the observed profile {(Jobs) for 
instrumental {oinst), thermal ((Tt/i), and fine struc- 
ture (cT/s) broadening. 



>th 



(1) 



We adopted an instrumental profile of FWHM of 
2.6 pixels (measured: 2.5-2.7 pi xels) for all ob- 



jects. We use the usual formula ( Lang 1980l eq 



2-243) to compute the thermal broadening, adopt- 
ing rest wavelengths of 6562. 83A and 5006. 85A 
for Ha and [O III] A5007, respectively, assuming 
no turbulent velocity, and adopting the observed 
electron temperature, when available (preferably 
from [O III] lines, but [N II] otherwise). If no elec- 
tron temperature was available, we used the mean 
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temperature for the other objects in each sam- 
ple (weak [O III] A5007 or hi gh He II 4686) For 
the fine structure broadening ( Meaburn1ll97Cll ). we 
adopted ct/s = 3.199 km/s for Ha and zero for 
[O III] A5007 (|Garcia-Diaz et ah Il2008al) . 

For real, spatially-resolved objects, the result- 
ing line width, AV, 



AV = 2.3556f7t^ 



(2) 



will be a luminosity-weighted, projected velocity 
width for the mass of the zone containing the emit- 
ting ion (O^"'' or H+) enclosed within the spec- 
trograph slit. We adopt half of this intrinsic line 
width in velocity units 



(3) 



as our measure of the kinematics for each object. 

We obtain angular diameters (see Table [T]) by 
collapsing the spectra along the wavelength axis 
to produce one-dimensional Ha spatial profiles, 
and then measured the diameter at 10% of the 

The uncer- 



peak intensity (jRicher et al. 
tainty in the diameters is less than half a pixel 
(O'.'S) for the sample with weak [O III]A5007 and 
no more than a full pixel (0'.'6) for the sample with 
strong He II 4686 (the Ha profiles contained at 
least 320,000 and 71,000 photons, respectively). 

3. Results and Discussion 

We separate the Bulge pla i ietary nebulae stud- 



ied here and by Riche r et al. I (|2008f ) into four evo- 



lutionary groups based upon the properties of the 
central star. The sample with weak [O III] A5007 
presented here has central stars whose temper- 
atures are sufficiently low that only low ioniza- 
tion ions exist and so should be the least evolved. 
Then follo w the younger [O III] A5007-bright ob- 



jects from [Richer et all (|2008i He II A6560 ab 



sent), with central stars that are hot enough to 
ionize O"'', but not 0^+ or He"*". The central stars 
in both groups should be on the horizontal por- 
tion of their post-AGB evolutionary track in the 
H-R diagram. The evolved [O III ] A5007-bright 



objects from [Richer et al~l pOOSl He II A6560 



present) and our sample with strong He II A4686 
have the hottest, most evolved central stars, with 
He^"*" present to differing degrees. While we expect 
some overlap between the two groups, the strong 



He II A4686 should be biased to later evolutionary 
stages ( §2.ip . These hotter central stars, partic- 
ularly those with strong He II A4686, may have 
extinguished nuclear reactions and be fading to- 
wards the white dwarf domain. 

Had we assigned our planetary nebulae exci- 
tatio n classes instead of defining our four groups 



,e.g ., AUer 1956; Gurzadvan iDopita fc Meatheringham 

[l990. : ,^Rcid fc Parker 2010|), there would be little 
difference in practice. Our weak [O III] A5007, 
young and evolved [O III] A5007-bright, and strong 
He II A4686 groups correspond to very low, low- 
to-medium, medium-to-high, and high excitation 
classes, respectively. Thus, our separation into 
evolutionary groups should be adequate for our 
purposes. 

3.1. Ha and [O III]A5007 line widths 

We present the relation between the line widths 
{AVq.z) in [O III]A5007 and Ha in Fig. H As 
has been found previously for bright Bulge plan- 
etary nebulae, th ere is a near equality be tween 



these fine widths (jRicher et al. Il2009l l201o[ ). The 



weak [O III] A5007 data set extends this relation to 
smaller line widths while the sample with strong 
He II A4686 is well-mixed with previous data. 

The tendency of finding systematically smaller 
line widths in [O III]A5007 at the smallest Ha line 
widths in Fig. El may be due to ionization strat- 
ification (e.g., IWilson I Il950f ). Hydrodynamical 
models clearly predict that the innermost parts 
of the ionized shell expand more slowly than the 
majority of the matter during the earliest phases 



of a p lanetary nebula's evolu tion (|Villaver et al. 



Perinotto et al. I 120041 ) since the central 
is not fast enough to create a hot 



20021 - 

star's wind 
bubble. 

The other feature in Fig. [6] is the kink near Ha 
line widths of about 33 km/s. The objects that 
lie at higher velocities than the kink all have hot, 
evolved central stars. Presumably, the kink results 
from a drop in the projected outfiow velocity in 
[O III]A5007, since that would seem energetically 
more feasible near the time when the central star 
is running out of nuclear energy. Might the kink 
be the result of the loss of pressure from the cen- 
tral star's wind after nuclea r reactions cease (col- 



lapse of the hot bubble; e.g..lGarcia-Segura et al 



20061 : iGarcfa-Diaz et al. I l2o'o8bh ? Ahernativelv. 
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the kink might also be produced via the pas- 
sage o£_a^JJe_^^jonization_^X)nt (Figs. 3 and 



5, Schonberner et al. 2005al ). Hydrodynam- 



ical models do not clea r ly predict this behav- 



ior (IVillaver et al. 1 12002': 'Per inotto et al. 112004 : 
Schonberner et al. | [2005a. 20^^~ More observa- 
tions and detailed modeling of individual objects 
will be required to understand the origin of the 
kink. 

3.2. Evolution of the Nebular Kinematics 

Fig. [7] presents the line width distributions for 
the four groups of Bulge planetary nebulae. The 
sample with weak [O III] A5007 clearly differen- 
tiates itself from the young [O III] A5007-bright 
objects, which is the evolutionary phase follow- 
ing it. On the other hand, there is no noticeable 
difference between the kinematics of the evolved 
[O III] A5007-bright objects and those with strong 
He II A4686. 

Statistical tests bear out these visual im- 
pressions. The no n-parametric U-test (e.g.. 



Wall fc Jenkins I I2003L §5.4.3) indicates that the 



probability of obtaining the line width distribu- 
tions for the objects with weak [O III] A5007 and 
the young [O III] A5007-bright objects from the 
same parent distribution is only 1.0 x 10~* when 
M3-13 is included (the highest line width in the 
weak [O III] A5007 sample) or 1.4 x 10"^ if it is 
excluded. Thus, the line width distributions for 
these two samples of objects are clearly different 
statistically. Not surprisingly, there is no statis- 
tical evidence for any difference in the line width 
distributions for the evolved [O III] A5007-bright 
objects and those with strong He II A4686. 

The line width distribution for the sample 
with weak [O III] A5007 (Fig. is very similar 
to the distribution of envelope expansion veloci- 
ties i n AGB stars ( Lewis l[l99ll : iRamstedt et al. 
2006f ) . This is exactly what hydrodynamical mod- 



els suggest for the earliest stages of evolution of 
the nebular shell before the central star's wind 



has created a hot bubble ( Villaver et al. I 2002 ; 



Perinotto et al. Il2004t iSchonberner et al. 112005a . 
20071) . Thus, this group of objects contains plan- 



etary nebulae whose nebular shells have not yet 
been significantly accelerated by the passage of the 



Richer et al. 



i oniza tion front and stellar winds 

(|2008[ ) lave already argued that the young and 
evolved samples of [O III] A5007-bright objects 



have the line width distributions expected from 
hydrodynamical models if they correspond, re- 
spectively, to the phase when a well-developed 
ionization front has swept through the nebular 
shell and the phase when the hot bubble is ac- 
tively accelerating the nebular shell. 

The line width distributions for the two groups 
of most evolved planetary nebulae are indistin- 
guishable. This implies that most of the nebular 
mass is not dramatically decelerated as the central 
star's luminosity decreases near the extinction of 
nuclear reactions, but continues expanding in a 
momentum-conserving mode. 

3.3. Nebular diameters 

Fig. [8] presents the distributions of angular di- 
ameters for the weak [O III] A5007 and strong 
He II A4686 sample s as well as those studied by 



Richer et al. I (|2008r ). Although the sample with 
weak [O III] A5007 are the least evolved, it does 
not have the smallest size distribution. Its distri- 
bution of angular diameters is statistically differ- 
ent from that for the young [O III] A5007-bright 
objects, with the U-test indicating a probability of 
only 6.1 X 10""* of drawing the two distributions 
from the same parent distribution. 

On the other hand, the distribution of diame- 
ters for the strong He II A4686 sample is clearly 
shifted to larger sizes compared to that for the 
evolved [O III] A5007-bright objects. The U-test 
indicates that the probability of drawing the two 
distributions from the same parent population is 
1.3 X 10~^. Assuming that the two groups are at 
the same distance, the difference in the distribu- 
tions of angular sizes also argues that the objects 
with He II A4686 are the most evolved. 

Hydrodynamical models might help explain 
these results. Models predict maximum H/3 lumi- 
nositi es when the [O III]A5007 emission is strong 
(e.g., Schonberner et al. 20Q7t ). Our requirement 
that the weak [O III]A5007 objects have high lu- 
minosity could cause us to preferentially select 
objects closer than the Bulge, which would make 
them appear larger. Also, models predict impor- 
tant structural changes during the early evolution, 
depending upon the relative impo rtance of the ion- 



izatio n front and hot bubble (e.g.. lPerinotto et al 



20041 ) that could affect the sizes we measure. 
Whatever the reason, a lack of a strong corre- 
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lation between nebular diameter and e volutionary 



indicator s have been ob tained before (jChu et al. 
[1984; Gu rzadvanlllQQTh . 



3.4. H/3 luminosities 

In Fig. [SJ we plot the line width as a func- 
tion of H/3 luminosity for the four planetary neb- 
ula samples. Alt hough the H/3 lumi nosities are 
not very accurate ( Richer et al. Ilioosh . there is a 
very clear progression as a function of evolution- 
ary state, from the weak [O III] A5007 sample, 
with the smallest line widths and high luminosi- 
ties, through the young and evolved [O III] A5007- 
bright objects to the strong He II A4686 sample, 
which has large line widths and the faintest lumi- 
nosities. 

The mixing of objects from the different sam- 
ples in Figs. [T] and [HlEl is undoubtedly due to 
the different effects that can affect the expan- 
sion of the nebular shell. Although the range in 
masses of the progenitor stars is likely relatively 
small, given their ages, the range in metallic- 
ity is substantially larger (e.g.. lSahu et al. Il200d : 



Zoccali et al. 1 120081 ). At lower metallicity the en- 



velope expa nsion velocities for AGB stars are also 
lower (e.g., Wood et al. 1992: Marshall et al 



2OO4I: 'Mattsson et al. "2008t IWachter et al. 1120081 
Grocnew cgcn ct al. 2009). On the other hand, 
hydrodynamical models find larger accelerations 
of the nebular shell at lower metall icity due to the 



higher electron temperature ( Sch onberner et al. 
2005bl ). The progenitor mass should affect the 
nebular kinematics primarily via the central 
star mass and its strong influence upon the 
wind output and evolutionary time scale (e.g., 
Villaver et al. ll2002t|Perinotto et al~ll2004l) . Con- 
sidering that the nebular shells are not likely to 
always be spherical, projection effects will further 
affect the measured line widths. Thus, the disper- 
sion seen in Figs. [T] and [Bjn] would seem plausible 
from natural causes. 

4. Conclusions 

We have obtained kinematic data for two sam- 
ples of planetary nebulae in the Milky Way bulge, 
selected so as to include objects very early and 
late in their evolution f i)2.ip . We measure line 
widths for Ha and [O III] A5007 in most cases. 
We combine these data sets with that studied by 



Richer et al. ( 20081 ). We define four evolutionary 
groups, based upon the temperature of the cen- 
tral star, which allow us to study the kinematics 
of the nebular shell from the earliest phases until 
the central star ceases nuclear burning. 

Generally, we find a near equality of the line 
widths for the Ha and [O III] A5007 lines in any 
given object. Ionization stratification likely ac- 
counts for the deviations: The [O III] A5007 line 
widths are systematically smaller than the Ha line 
widths for the smallest and largest Ha line widths, 
corresponding to the earliest and latest evolution- 
ary phases, respectively. 

We see clear evolution of the kinematics of the 
nebular shell. The least evolved objects, our plan- 
etary nebulae with weak [O III] A5007, have cool 
central stars and the nebular envelopes have a 
line width distribution similar to that of the enve- 
lope expansion velocities of AGB stars, indicating 
that the ionization front has not yet been able to 
substantially ac celerate the nebular shell. In sub- 



sequent phases (jRicher et al. 1120081 ). the nebular 



shell is first accelerated by the passage of an ion- 
ization front and then further accelerated once the 
central star's wind produces a hot bubble. The 
line width distributions for the planetary nebu- 
lae in these three evolutionary phases are statisti- 
cally distinct. The most evolved objects, with high 
He II A4686 ratios, have a similar line width dis- 
tribution to evolved [O III] A5007-bright objects, 
suggesting that no further acceleration occurs as 
the central stars reach their highest temperatures, 
their nuclear reactions cease, and their winds de- 
cline. 

This kinematic evolution of the nebular shell 
has lon g been predicted by hydrody namical mod^ 
els I 



e.g 



20m 



Kahn fc West 111985: Marten & Schonberner 



199ll : lMellemalll994HVillaver et al.'1l2002HPerinotto et al 



Richer et al 



Our resul ts, together with those of 
( 2008[ ). based upon a large sam- 



ple of Galactic planetary nebulae from a single 
stellar population, clearly confirm these predic- 
tions. At least until the point at which nuclear 
reactions cease in the central stars, their ionizing 
fiuxes and winds continuously accelerate the neb- 
ular envelopes. What happens thereafter is not 
yet completely clear, and would require samples 
of planetary nebulae chosen specifically to contain 
hot central stars of low luminosity. 
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Table 1 

Bulge Planetary Nebula Sample 



Objocl. 


PN G 


sample 


Run 


FWHMfHaV' 


AV„ -,(Ha) 


FWHM(5007)^' 


AV„ ,r,(5007) 


H(v Diameter'' 














(A) 


(klll/s) 




(A) 




^kln/^j 






H 1-24 


004.6-1-06.0 


weak 5007 


2008jun 





0700 _|_ n C\C\dA 


11.22 ± 0.10 


0. 


ftOQQ _|_ n 


,0081 


17.54 ± 0.24 


7.3 


34.84 


H 1-34 


005.5-1-02.7 


weak 5007 


2008jun 





SI fiS -1- n 0079 


12.31 ± 0.16 


0. 


KAOO -U n 


,0077 


14.88 ± 0.23 


5.0 


34.60 


H 1-39 


356.5-03.9 


weak 5007 


2008jun 





7707 _|_ n 001 Q 


11.31 ± 0.04 


0. 




,0030 


7.77 ± 0.09 


5.9 


34.80 


H 1-43 


357.1-04.7 


weak 5007 


2008jun 





70Qn -1- 00*^0 


10.82 ± 0.07 










5.8 


34.45 


H 1-44 


358.9-03.7 


weak 5007 


2008jun 





«-! QQ _|_ n 0099 


5.37 ± 0.05 


Q 




0028 


6.19 ± 0.08 


6.3 


34.10 


H 1-55 


001.7-04.4 


weak 5007 


2008jun 


Q 


f^Qcn _!_ n nOI 4 

.\JO\Ji} ZL. \J .\>\J L'-t 


6.30 ± 0.03 


Q 


2939 ± 


0047 


6.18 ± 0.14 


5.7 


34.38 


H 2-25 


004.8+02.0 


weak 5007 


2008jun 


Q 


.VJfiiO ZL. \J .\>\J ^0 


2.37 ± 0.05 


Q 


2747 ± 


0087 


4.79 ± 0.26 


7.1 


34.54 


H 2-29 


357.6-03.3 


weak 5007 


2008jun 


Q 


QiA'i'i X- 00R9 


15.28 ± 0.14 










12.7 


33.49 


H 2-48 


011.3-09.4 


weak 5007 


2008jun 


Q 


^^aA -V- 001 9 


9.45 ± 0.03 





Q9Q0 _|_ n 


0015 


7.42 ± 0.04 


5.9 


35.26 


He 2-260 


008.2-1-06.8 


weak 5007 


2008jun 





7001 _|_ n 001 2 


6.48 ± 0.03 


0. 




,0035 


5.91 ± 0.10 


7.7 


34.39 


M 1-26 


358.9-00.7 


weak 5007 


2008jun 





70E0 -L n 009^ 


10.13 ± 0.05 


0. 


.3477 ± 0. 


,0031 


8.23 ± 0.09 


7.3 


36.07 


M 1-27 


356.5-02.3 


weak 5007 


2008jun 





7970 -L n 001 S 


10.43 ± 0.04 










8.9 


35.54 


M 1-30 


355.9-04.2 


weak 5007 


2008jun 





0701 -1- 009Q 


14.10 ± 0.07 


0. 


554Q + 


,0058 


15.32 ± 0.17 


6.5 


34.97 


M 1-44 


004.9-04.9 


weak 5007 


2008jun 





p.'7Af\ -i- n ooso 


8.13 ± 0.07 










7.6 


34.75 


M 1-45 


012.6-02.6 


weak 5007 


2008jun 


Q 


77CQ J_ n n098 

. 1 / VJO ZL. \J .\>\J ^0 


12.15 ± 0.06 










5.6 


35.15 


M 2-10 


354.2-1-04.3 


weak 5007 


2008jun 


Q 


fiOT c _|_ n nni n 

.Uf -LU ZL. \J .\j\J i-\J 


9.00 ± 0.02 


Q 




0024 


10.21 ± 0.07 


7.5 


34.90 


M 2-12 


359.8+05.6 


weak 5007 


2008jun 


Q 


f;-] Qc _|_ n noi 1 


4.81 ± 0.03 


Q 


2955 Q 


0447 


6.13 ± 1.34 


7.1 


34.56 


M 2-14 


003.6+03.1 


weak 5007 


2008jun 


Q 


77n'i -1- n 

. 1 (UO ^ U.UUOO 


11.30 ± 0.09 


Q 


AA^fi 4- n 

,'±'±^K) HZ U. 


0048 


11.63 ± 0.14 


5.3 


34.47 


M 2-19 


000.2-01.9 


weak 5007 


2008jun 





7C04 _|_ n 0099 


10.59 ± 0.05 


0. 


. ooou m o. 


,0036 


9.73 ± 0.11 


8.8 


34.75 


M 2-7 


353.7+06.3 


weak 5007 


2008jun 





0400 -1- n oo9fi 


14.13 ± 0.06 


0. 


. uoo 1 m o. 


,0035 


14.61 ± 0.10 


9.9 


34.15 


M 3-13 


005.2+04.2 


weak 5007 


2008jun 




oo7fi _|_ n 09S9 


26.20 ± 0.64 


0. 


. ^ i \jo m o. 


,0556 


28.38 ± 1.67 


6.2 


35.24 


M 3-17 


359.3-03.1 


weak 5007 


2008jun 





001 e _|_ n 0090 


13.47 ± 0.05 


0. 


.4159 ± 0. 


,0051 


10.70 ± 0.15 


5.7 


34.90 


SwSt 1 


001.5-06.7 


weak 5007 


2008jun 


Q 


Ck077 -\- n 0097 


16.43 ± 0.06 


Q 


nqi q -L n 


0033 


16.54 ± 0.10 


4.6 


35.03 


Th 3-16 


357.5-03.1 


weak 5007 


2n08j un 


0. 


fiCQo _L n noi 7 


7.10 ± 0.04 










6.3 


34.05 


Al 1 


006.8-08.6 


strong 4686 


2008sep 


1 


9709 _|_ n nO'^fi 


24.04 ± 0.17 


Q 




0162 


21.73 ± 0.49 


13.0 


34.12 


Al 2-E 


359.3+03.6 


strong 4686 


2008sep 


1 




33.31 ± 0.74 




0817 i 


0283 


31.66 ± 0.85 


8.6 


33.19 


Al 2-H 


357.2+01.4 


strong 4686 


2008sep 


2 




16.07 ± 0.19 


Q 


7402 + 


0180 


21.08 ± 0.54 


9.3 




Al 2-1 


359.5+02.6 


strong 4686 


2008sep 




iQcn _|_ n OOQQ 


21.41 ± 0.30 


0. 




,0210 


26.25 ± 0.63 


6.6 




H 2-44 


005.5-04.0 


strong 4686 


2008sep 


\ 


AQ-l 9-1-0 09fifi 


29.51 ± 0.80 


1 


.1545 ± 0. 


,0395 


33.89 ± 1.18 


12.4 


34.16 


KFL 02 




strong 4686 


2008sep 


1. 


.2658 ± 0.0102 


23.70 ± 0.31 


0. 


,9265 ± 


,0169 


26.89 ± 0.51 


8.6 


31.49 


KFL 09 




strong 4686 


2008sep 


1. 


.3491 ± 0.0240 


25.99 ± 0.72 


1, 


,0010 ± 0, 


,0259 


29.18 ± 0.78 


12.0 


33.82 


KFL 16 


UU0.D-U4. ( 


strong 4686 


o n n o r-i ^ 
zUUoSCp 


1. 


.5910 ± 0.0254 


OZ.OD ± U. (o 


1, 


,2402 ± 0, 


,0450 


06.OI ± 1.00 


13.9 


33.38 


M 2-38 


005.7-05.3 


strong 4686 


2008scp 


1. 


.3917 ± 0.0185 


27.95 ± 0.55 


0, 


,9841 ± 0, 


,0315 


28.71 ± 0.94 


14.4 


33.50 


M 3-22 


000.7-03.7 


strong 4686 


2008scp 


1. 


.2684 ± 0.0125 


23.70 ± 0.37 


0, 


,9609 ± 0, 


,0250 


27.94 ± 0.75 


9.3 


34.21 


M 3-23 


000.9-04.8 


strong 4686 


2008scp 


1. 


.4720 ± 0.0291 


29.27 ± 0.87 


1, 


,0804 ± 0, 


,0538 


31.62 ± 1.61 


14.0 


34.69 


Pe 1-12 


004.0-05.8 


strong 4686 


2008sep 


1, 


.5601 ± 0.0428 


31.84 ± 1.28 


1, 


,2239 ± 0, 


,0589 


36.03 ± 1.76 


11.9 


33.38 


Pe 1-13 


010.7-06.7 


strong 4686 


2008sep 


1, 


.0811 ± 0.0059 


18.42 ±0.18 










9.4 


33.69 


Pe 2-13 


006.4-04.6 


strong 4686 


2008sep 


1, 


.3648 ± 0.0129 


26.42 ± 0.39 










9.8 


34.08 


Sa 2-230 


010.7+07.4 


strong 4686 


2008sep 


1, 


.7239 ± 0.0576 


35.74 ± 1.73 


1, 


,1654 ± 0, 


,0490 


34.23 ± 1.47 


15.1 


33.63 


SB 15 


009.3-06.5 


strong 4686 


2008sep 


1, 


.2971 ± 0.0067 


23.40 ± 0.20 


0, 


,9568 ± 0, 


,0116 


27.76 ± 0.35 


8.5 


32.76 


SB 37 


352.6-04.9 


strong 4686 


2008scp 


1. 


.5842 ± 0.0175 


32.18 ± 0.52 


1, 


,2002 ± 0, 


,0297 


35.29 ± 0.89 


10.0 


33.92 


SB 38 


352.7-08.4 


strong 4686 


2008scp 


1. 


.4677 ± 0.0154 


28.59 ± 0.46 


1, 


,1736 ± 0, 


,0255 


34.43 ± 0.76 


14.0 


32.18 


SB 55 


359.4-08.5 


strong 4686 


2008sep 


1, 


.4107 + 0.0356 


25.56 ± 1.07 


1, 


,0770 ± 0, 


,0490 


31.41 ± 1.47 


14.9 


32.88 


ShWi 2-1 


001.4-03.4 


strong 4686 


2008sep 


1, 


.0534 ± 0.0250 


17.45 ±0.75 










11.8 


33.60 



Table 1 — Continued 



Object 


PN G 


sample 


Run 


FWHM(Ha)* 

(A) 


AVb.5(Ha) 
(km/s) 


FWHM(5007)'' 

(A) 


AVo.5(5007) 
(km/s) 


Ha Diameter^ 

10% Imax 


(erg/s) 


Th 3-26 
H 2-13 
H 2-43 


358.8+03.0 
357.2+02.0 
003.4-04.8 


strong 4686 
6560 present 
6560 absennt 


2008sep 
2008jun 
2008jun 


1.0768 ± 0.0090 
0.9376 ± 0.0023 
1.9363 ± 0.0185 


20.97 + 0.27 
16.72 + 0.05 
39.20 + 0.42 


0.6937 + 0.0103 
0.5785 + 0.0045 
1.7300 + 0.0658 


19.78 + 0.31 
16.12 + 0.13 
51.28+ 1.97 


9.7 

6.0 
4.6 


34.05 
34.47 
34.88 



^This is the observed line width, uncorrected for instrumental, thermal, or fine structure broadening. 

^The diameter was measured by collapsing the Ha spectra into one-dimensional spatial profiles and measuring the diameter at 10% of the maximum 
intensity. 
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Fig. 1. — The selection criteria for the two sam- 
ples of Bulge planetary nebulae presented here 
were designed to select planetary nebulae that 
are (a) less evolved a nd (b) more evolved than 
the sample studied bv iRicher et al. I (|2008l ). The 



area within the dotted box is approximately the 
area occupied by bright extragalactic planetary 



tion 


( 


Jacobv & Ciardullo 19991 Richer et al. 


1999 


• 


Walsh et al. 


19991 Dudziak et al. 


200C 


■ 


Roth et al. 1 12004 


: Mendez et al. 1 


2005; 


Ziilstraa et al. l2006l 


Goncalves et al. 


2007; 


Richer & McCall 112008). 
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Fig. 7. — There is a very clear increase in the 
line width, progressing from planetary nebulae 
with weak [O III]A5007 to the young and evolved 
[O III]A5007-bright objects, clearly demonstrating 
that the nebular shells are accelerated throughout 
these phases. From the perspective of line widths, 
these three phases are statistically distinct. Plan- 
etary nebulae in the most evolved phase, with 
strong He II 4686 line intensities, have line widths 
very similar to the evolved [O III] A5007- bright 
planetary nebulae. M3-13 is the "anomalous" ob- 
ject with weak [O III]A5007 and a line width of 
26 kr n/s (see Fig. [3| . We correct the line widths 



i^ig. m. w 

~ari (|2008 



from lRicher et al. I (|2008l ) for thermal broadening 
using the observed electron temperature, instead 
of a constant value of 10'* K. As a result, our his- 
tograms for those data differ. 
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Fig. 2. - We present the line profiles in [O III]A5007 and Ha for the sample of planetary nebulae with 
weak [O III]A5007 lines. The wavelength interval plotted is always that indicated on the bottom panels. 
The number beside the name in the Ha panels is the angular diameter of the object at 10% of maximum 
intensity. 
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5000 5003.5 5005 5007-5 5010 5012.5 5015 6555 6560 6565 6570 6575 

lavslenglh (aiiESlroms) lavelenglh (alisalroms) 



Fig. 3. — As in Fig. we present the line profiles for the sample of planetary nebulae with weak [O III]A5007 
lines. Note the very wide wings on the profiles for M 3-13. 
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Fig. 4. — As in Fig. [51 we present the line profiles for the sample of planetary nebulae with strong He II 
4686 lines. 
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Fig. 5. — As in Fig. [51 we present the line profiles for the sample of planetary nebulae with strong He II 
4686 lines. 
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Fig. 6. — The intrinsic line width in [O III]A5007 is plott ed as a function of th e intrinsic hne width in Ha for 



The intrinsic hne widths are 



the two samples presented here as well as the data from [Richer et al. 
corrected for the instrumental resolution, fine structure broadening, and thermal broadening. The solid line 
indicates the locus of identical line widths in the two emission lines. The weak [O III] A5007 sample extends 



the results of [Richer et al. I (|2008f ) to lower line widths while the strong He II A4686 sample is well mixed 
with their data. The error bars appear to fill many of the symbols for the objects in the samples with weak 
[O III] A5007 and He II A6560 absent. 
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Fig. 9. — The line widths are plotted as a function of the H/3 luminosity for the four samples of Bulge 
planetary nebulae. A variety of effects are responsible for mixing the objects from different evolutionary 
stages. 
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Fig. 8. — There is a clear increase in the diam- 
eters for the planetary nebulae in the young and 
evolved [O III]A5007-bright and strong He II 4686 
groups. The least evolved planetary nebulae, with 
with weak [O III]A5007, break this evolutionary 
sequence. 
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